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Magnetic State Modification Induced by Superconducting Response in 
Ferromagnet /Superconductor Hybrids 
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Magnetization measurements in superconductor/ferromagnet Nb/Co superlattices show a com- 
plex behavior as a function of temperature, applied field and sample history. In base to a simple 
model it is shown that this behavior is due to an interplay between the superconductor magne- 
tization temperature dependence, the ferromagnet magnetization time dependence, and the stray 
fields of both materials. It is also shown that the magnetic state of the Co layers is modified by 
the A^6 superconducting response, implying that the problem of a superconductor/ferromagnetic 
heterogeneous sample has to be solved in a self-consistent manner. 



The interaction between a superconductor, SC, and a 
ferromagnet, FM, in close contact at an interface, as in a 
superlattice, has attracted attention in the last years due 
to the possibility of fabricating SC/FM hybrid devices. [H 
These engineered materials originate the appearance of 
interesting physical phenomena due to the different scales 
and mechanisms of interaction, like SC pair breaking ef- 
fects related to exchange interaction at the interface, [1] 
or electromagnetic interaction with the stray fields of the 
FM both at the mesoscopic and macroscopic level. 0, Q 

Most of the research has focused on the ways in 
which the FM affects the SC response. For example, 
in the Domain Wall Superconductivity effect, observed 
in SC/FM bi-layers,[l, 13| superconductivity nucleates 
in those places where the perpendicular component of 
the inhomogeneous FM domain structure's stray field 
is close to zero. In the Spin Switch effect, observed in 
FM/SC/FM trilayers,[ll the Cooper pair, due to its fi- 
nite size, experiments different average values of the ex- 
change field when the FM layers are ferro- or antiferro- 
magnetically oriented. As a consequence, the SC order 
parameter is more depressed when the FM layers are 
ferro-magnetically oriented. This allows the control of 
the SC order parameter value through an external macro- 
scopic parameter. 

In contrast, very little work has been done in explor- 
ing in which way the SC affects the magnetic state of 
the FM layer. Q Recently, we have shown the impor- 
tance of the FM stray fields in the overall magnetic re- 
sponse of Nb/Co superlattices. In that work we also 
hinted to the possibility that the SC response may mod- 
ify the magnetic state of the FM layers. In this letter we 
show that, indeed, the SC response modifies the magnetic 
state of the FM layers. The system global electromag- 
netic response is determined by an interplay between the 
SC magnetization temperature dependence and the FM 
magnetization time evolution. 

We present data on the temperature, T, dependence 
of the magnetic flux expulsion, A(f), directly propor- 
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FIG. 1: Main panel: Magnetization, M, as a function of ap- 
plied field. Ha, for the [A''6(44 nm)/Co( 10 nm)]xl9 superlat- 
tice at T = 7iC > Tcs = 6.2 K. Panels (a) through (d) 
show the temperature, T, dependence of the superconducting 
magnetic fiux response, A(j>, in units of the superconducting 
flux quantum, (/lo, at different applied fields and ferromagnetic 
layers initial state, as indicated by gray dots and connecting 
arrows in the main panel. Arrows in panels indicate the direc- 
tion of the T sweeps. Panel (a): Ha = -22 Oe, -l-FC. Panel 
(b): Ha = 75 Oe, -FC. Panel (c): Ha = 38 Oe, -FC. Panel 
(d): Ha = 22 Oe, -FC. 



tional to the SC magnetization, for Nb/Co superlattices. 
Data is presented for a [Nb{4Anm)/Co{10n'm)]xl9 and 
a [A^6(44nm)/Co(7.5nm)]xI9 superlattice. The results 
on these samples are representative of the measurements 
we performed on a collection of Nb/Co FM/SC super- 
lattices. Since our experimental setup measures mag- 
netic flux variations,!^ all Ac/) values are measured with 
respect to the first data point, always above the super- 
conducting critical temperature, Tcs- Sample prepara- 
tion, characterization method and measurement details 
are described in reference [HI. For all flux expulsion data. 
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the applied field, Ha, is parallel to the sample surface. 
In the normal state, both superlattices present FM be- 
havior with Curie temperatures, Tq, above 300 Flux 
expulsion in the SC state was measured as a function of T 
in field cooling experiments, for two different Co layer's 
initial FM states, -FC, with the Co layers initially sat- 
urated in the negative Ha direction, and -|-FC, with the 
Co layers initially saturated in the positive Ha direction. 
A detailed explanation of these measurement protocols 
is also included in reference 0. 

Figure [1] is a composite that summarizes the experi- 
mental results. In the main panel we show the FM hys- 
teresis loop of the Co layers at T = 7 K, close but above 
Tcs of the Nb layers. In the superimposed panels, we 
show the T dependence of A(f) for -|-FC and -FC mea- 
surements. For each initial state, several T cycles were 
measured, each cycle sweeping T down from 7 K to 5.5 K 
and up to 7 K again. The solid dots in the hysteresis 
curve connected with arrows to the panels indicate the 
initial magnetic state for each experiment. 

The data for the first T down-sweep in each panel 
shows the behavior already discussed in our previous 
work.0. As discussed there, the SC response is pro- 
portional to the effective field, H^ff, originated by the 
superposition of the applied field. Ha, and the Co lay- 
ers' stray field, Hg. For the -FC measurements (lower 
branch of the Co hysteresis loop), at low Ha and neg- 
ative Co magnetization, the SC layers sense a positive 
Heff due to the Co's Hs, see panel (d). At higher Ha, 
the Co magnetization becomes positive, Hg becomes neg- 
ative and larger than Ha, and the SC senses a negative 
Hgff, see panel (b). Panel (c) shows an intermediate 
case, where the magnetization is already reversed, but 
Hs is smaller than Ha and the SC still senses a positive 
Heff- Panel (a), -l-FC initial state (upper branch of the 
Co hysteresis loop) is the mirror experiment from panel 
(d). 

The novel feature observed in these data is present in 
the dependence of the Co magnetization in the normal 
state with the number of T sweeps, i.e. cycles. This de- 
pendence is observed as a non-repeatability of the normal 
state magnetization value after a cycle is completed. This 
behavior is not due to an experimental artifact related to 
an instrumental drift, since this instrumental drift has 
been substracted from the data. A systematic behavior 
is observed in spite of the seemingly complex dependence. 
The direction of the variation follows the sign of the ap- 
plied field, and is independent of the Co magnetization 
direction, i.e. the stray field, compare data in panel (b) 
and (d), for example. The difference between the first 
and second normal state A(/) values in panel (d) of 360 
superconducting flux quantums is equivalent to a change 
of 1.1 emu cm^^ in the Co layer's magnetization, which 
shows that this effect is small but not negligible. 

In order to understand this behavior, we have con- 
structed a simple "toy model" to qualitatively simulate 




FIG. 2: Schematics of the "toy model" behavior. Panel (a): 
at T > Tcs the normal A''6 ellipsoid experiences an effective 
field due to the applied field (straight lines) and the ferromag- 
netic Co ellipsoid stray field (dipole like lines). Panel (b): As 
T is reduced below Tcs the magnetic flux expulsion from the 
Nb ellipsoid modifies the effective field over the Co ellipsoid. 



the experimental data. Although the model is very sim- 
ple, a careful consideration of its hypothesis should be 
made to fully understand the implications of the results. 
The first requirement is that an electromagnetic stray- 
field mediated interaction should exist between the FM 
and the SC components. This is not achievable if the 
materials are modelled as nearly infinite slabs parallel to 
the applied field, since the stray field of this geometry 
is negligible. Consequently, both materials are modelled 
as ellipsoids with one of the principal axis parallel to the 
field. The "toy model" sample consists, then, of a FM 
and a SC ellipsoids, located side by side. Figure [2] de- 
picts the main ideas of the model. Panel (a) shows the 
situation at T > Tcs where the Heff sensed by the SC 
ellipsoid is composed by Ha (straight lines) and the FM 
Hs (dipolar lines arising from the FM ellipsoid) . When T 
is reduced below Tcs, depicted in panel (b), the flux 
expulsion from the SC ellipsoid modifies the Hef / sensed 
by the FM ellipsoid and consequently its magnetization. 
The solution of the problem has now to be found in a 
self-consistent way. 

The ellipsoid shape or eccentricity, e, was selected as 
to maximize the stray field effects. That an optimum 
value exists is clear from the fact that in the e — > cxd 
"needle" limit, the stray fields approach zero due to the 
negligible demagnetizing effects, and that in the e — > 
"disk" limit, the stray fields also approach zero since the 
ellipsoid is being magnetized along the shape anisotropy 
"hard axis" . The optimal e value actually depends on the 
material's magnetization, but since it is weakly depen- 
dent on it, a value of 10 was found to maximize the stray 
field effects in nearly all the T-Ha range. Also, since an 
exact three dimensional spatial solution of this electro- 
magnetic problem is beyond the scope of this work, and 
would only obscure the results of the model, the spatial 
dependence of the stray fields is neglected, and Hs due to 
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each ellipsoid is evaluated only at the center of the other 
ellipsoid. 

The second ingredient in the model is a "time" de- 
pendence. This dependence cannot be ascribed to the 
superconducting material since we have shown that no 
vortices are present in the T — Ha range of these exper- 
iments. Consequently, it must be arising from the creep 
in the FM material. Following this idea, the magnetiza- 
tion of the SC Nb ellipsoid is modelled by a T dependent, 
time independent Meissner state. As a further simplifi- 
cation of the model the T dependence is forced to follow 
that of a parallel slab with a two fluid T law. Q On the 
other hand, the FM Co ellipsoid magnetization does not 
present a T dependence since its Tc is much higher than 
the measurement range. It only shows a time dependence 
which must be numerically simulated, as described in the 
next paragraph. 

To simulate the T sweeps at constant Ha, the self- 
consistent equilibrium state of the magnetized ellipsoids 
is solved at a given T. After this, the magnetization 
change for the Co ellipsoid, is reduced by a given per- 
centage, and the magnetization of the Nb ellipsoid is 
recalculated for this, now fixed, value of the Co ellip- 
soid's magnetization, i.e. stray field. This algorithm re- 
sults in an effective exponential time dependence for the 
Co magnetization. The sample's magnetization, Mt, 
is defined as the total magnetic moment divided by the 
total sample volume. In order to compare the results 
to the experiments, the simulation data is presented as 
AMt = Mt — Mo, where Mo is the value for the first 
simulated point, always at T > Tcs- 

Panels (a) and (b) in figure [3] show the prediction of 
the model for situations similar to panels (c) and (a) in 
figure [H i.e. opposite direction of Ha and same value of 
Mo- It is clear that the principal features of the exper- 
imental data are qualitatively reproduced. First, there 
is a dependence of the normal state magnetization with 
the number of cycles. This dependence follows the sign 
of the applied field and is not correlated to the mag- 
netization direction. Second, there is an irreversibility 
between cooling-down and warming-up sweeps. Third, 
a non-monotonic T dependence is observed for cooling- 
down sweeps. 

An interesting feature not actually observable in the 
data in figdj but presented in reference d is a non- 
monotonic T dependence that develops for applied fields 
near the coercive field of the Co layers. The main panel in 
figure [4] shows a comparison between experimental and 
simulated data, where the simulation parameters have 
been selected as to maximize this non-monotonic T de- 
pendence. The origin of this behavior becomes clear 
when examining separately the Nb and Co magnetiza- 
tion response in the simulated data. The inset shows the 
SC ellipsoid magnetization, Msc, the FM ellipsoid mag- 
netization, MpM, and the sample's magnetization, Mt, 
as a function of simulated data point number, i.e. "time" , 
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FIG. 3: Model prediction for the temperature, T, dependence 
of the sample magnetization change, AMt = Mt{T) — Mo, 
where Mo is the magnetization for the initial simulated data 
point at T > Tcs- Panels (a) and (b) show the results 
for two sets of parameters qualitatively equivalent to the ex- 
perimental data in panels (c) and (a) of figure [1] Panel 
(a): Ha = 22.85 Oe, Mo = 5.12 em«cm~^ Panel (b): 
Ha = -17.15 Oe, Mo = 5.12 emucm"^ 



while T is swept down from above Tcs- The T sweep is 
linear with this "simulated time" . The time dependence 
of AIsc is that arising from the T sweep, given that the 
Meissner state does not present an intrinsic time depen- 
dence. The time dependence of MpM, on the other hand, 
has a twofold origin. First, the fiux expulsion in the SC 
originates an increase of local magnetic field in the FM 
material, as schematized in panel (b) of figure[2] Second, 
the MpM presents an intrinsic time dependence in its re- 
sponse to the magnetic field changes. In this light, the 
origin of the non-monotonic T dependence becomes clear. 
As T is swept down from above Tcs, the T dependence 
of the SC ellipsoid magnetization produces a flux expul- 
sion in the sample. This originates a field increase in the 
FM material increasing its MpM- At lower temperatures, 
the T dependence of the SC material is relatively weak, 
and the time dependence of the FM material emerges as 
a "paramagnetic" like signal, resembling a paramagnetic 
Meissner effect, jl, 0, [1]. 

The results and the toy model presented here clarify 
the response of SC/FM hybrid structures and, at the 
same time, raise an interesting question. We have demon- 
strated that the electrodynamic response of these hybrid 
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FIG. 4: Magnetic flux expulsion, for a 

[Ar6(44 nm)/Co(7.5 nm)]xl9 superlattice with Tcs = 5.9 K, 
and model prediction, AMr, for Ha = 15 Oe and 

Mo = Qemucm^"*. Inset shows the dependence on 
simulated data point number, i.e. "time", for the supercon- 
ducting ellipsoid magnetization, Msc, the ferromagnetic 
ellipsoid magnetization, Mfm, and the sample's total 
magnetization, Mt- 

systems involves a combination of two separate phenom- 
ena. In the first place, the diamagnctic response of the 
SC layers expels the magnetic flux into the FM layers. 
As a consequence, the FM material responds with a time 
dependence, clearly in the direction of the applied field. 
Both materials affect each other with their respective 
stray fields. In this process, the magnetic domain struc- 
ture of the FM seems to play an important role, since 
the stray fields of an infinite slab are negligible. Clearly, 
in order to observe stray field effects, a non slab geome- 
try has to be present in the samples. In this picture, an 



interesting point arises. Given that the response of the 
hybrid material is affected, and in some T and H range, 
dominated, by the intrinsic time dependence, the effects 
described here may be important if the device operation 
is based on magnetization changes and designed to work 
at frequencies similar to the creep of the FM. 

In summary, we have demonstrated that the electro- 
dynamic response of SC/FM hybrid materials is deter- 
mined by an interplay between the temperature depen- 
dence of the SC magnetization, the time dependence of 
the FM magnetization, and the effective interaction be- 
tween them mediated by the stray fields. 
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